Domain walls in supersymmetric QCD: from weak to strong coupling 
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We consider domain walls that appear in supersymmetric QCD with Nf < Nc massive flavours. 
In particular, for 2Nf < Nc we explicitly construct the domain walls that interpolate between vacua 
labeled by i and i + Nf. We show that these solutions are Bogomol'nyi-Prasad-Sommerfield (BPS) 
saturated for any value of the mass of the matter fields. This fact allows us to evaluate the large mass 
limit of these domain walls. We comment on the relevance of these solutions for supersymmetric 
gluodynamics. 
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In recent times, a lot of attention has been drawn to 
the existence of exact solutions in supersymmetric gauge 
theories that are in the strong coupling regime. One of 
the more relevant issues is that of domain walls in SU(N) 
supersymmetric gluodynamics, the theory of gluons and 
gluinos. Those arise because this theory has an axial 
U(l) symmetry broken by the anomaly to a discrete Z2N 
chiral symmetry. Due to non-perturbative effects gluino 
condensates ( (AA) ) form, breaking the symmetry further 
down to Z2 ■ This leaves us with a set of N different vacua 
labelled by 



(TrAA) = A3^27rik/N ^ 0, 1, ...,N - 1 



(1) 



where A is the condensation scale, and, as indicated 
above, a set of domain walls interpolating between them 
(in Ref. it was pointed out the existence of a chirally 
symmetric vacuum where the gaugino condensate van- 
ishes. Here we will only consider domain walls involving 
chirally asymmetric vacua). If we asume that they are 
BPS saturated, the energy density of these walls is ex- 
actly calculable and given by 



e = 



TV 
8^ 



|(TrAA>oo - (TrAA)_c 



(2) 



In fact, it has been suggested in Ref. |^ that, in the large 
N limit, these domain walls are BPS states. On the other 
hand, these solutions preserving half of the supersymme- 
try would play an important role in the D-brane descrip- 
tion of = 1 SQCD 0. However, we want to stress that 
whether or not these configurations were BPS saturated 
was, up to now, still an open question 

A useful way of gaining intuition on pure gluodynamics 
is by adding matter fields and analyzing the limit where 
these extra fields become very heavy. These new fields 
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are usually taken to be pairs of chiral superfields trans- 
forming as {N, N) under the color group. In the strong 
coupling regime, squark condensates will thefore form. 
These models, for the case of (A^ — 1) flavours, were con- 
sidered in Refs |8|-[lT| , where the analysis of the vacuum 
structure led the authors to conclude that the existence 
of BPS saturated domain walls was restricted to values 
for the mass m of the squark fields below a certain critical 
one. This jeopardized the idea of recovering pure gluo- 
dynamics by taking the limit m — > 00, which is precisely 
the issue we want to revise here. 

In order to do that let us consider supersymmetric 
QCD with SU{Nc) gauge group and Nf couples of chi- 
ral superfields {Qi,Qi) transforming as Nc,Nc- Non- 
perturbative effects become relevant at the scale A, where 
condensates form. The gaugino and squark colorless con- 
densates are described by the following composite fields 



S 



327r^ 



(3) 



M'^^Q'Q, i,j^l,2,...,N 



f 



where is the composite chiral superfield whose low- 
est component is AA. In this regime, the relevant degrees 
of freedom are described by a Wess-Zumino model, as 
shown in Ref. [n2| . Its effective Langrangian is given by 



(fe w + h.c. 



(4) 



where JC is the Kahler potential and W is the superpo- 
tential 



2 ^^c-iV/detM 

3 A3^<=-^/e^--^/ 
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Tr(TOAf) 



(5) 
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with mi the mass matrix for the matter superfields. 



We 



will work in the flavour basis where this matrix is diago- 
nal, and will use the notation m-^ = Slrrij. We will ana- 
lyze the simple case where the corresponding eigenvalues 
are real. This superpotential has Nc extrema labeled by 
the different phases of the gaugino condensate. At the 
minimum we have the gaugino condensate fixed to 



— detm 
4/ 



(6) 



where, and from now on, we set A = 1. The matter fields 
are aligned with respect to the former and given by 



ML 
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Finally, the superpotential at the minimum is propor- 
tional to the gaugino condensate 



(8) 



We want to study domain wall configurations that in- 
terpolate between the different minima. Here a technical 
problem appears: the superpotential has several branches 
associated with its logarithmic piece p^ . In the pure 
SUSY gluodynamics limit described by Veneziano and 
Yankielowicz Jl2[ this is a severe problem, since any con- 
figuration going from one vacuum to another has to cross 
this branch. This is not necessarily the case when we in- 
clude other fields, given that the variation in the phase of 
the gaugino condensate can be partially compensated by 
these new fields. In this case, this will be done by matter 
fields. 

Let (S", M)a be a particular vacuum. We can continu- 
ously deform it into another vacuum, {S,M)i,. For this 
path in the configuration space, we define S, Wi such that 



S\b = e''S\a 



(9) 



Since, as mentioned above, matter at the minimum has 
to be aligned with respect to the gaugino condensate, 
Wi must be some integer numbers. On the other hand, 
one necessary condition to avoid crossing the logarithmic 
branch along this general path is 



(Nc - Nf )6 + J2iS + 2nw,) = , 



(10) 



and then 6 = 27r-|^, where k is the integer given by 

If we assume that there is a BPS domain wall connect- 
ing these two vacua from z = — oo to z = -foo, it will be 
described by the following differential equations 



(11) 



KMMdM = e-^ 



aw 

dM] 



c2 tr 

where K,^^ = is the induced metric from the Kahler 
potential /C, and 7 is given by 



1 , J. , kn TT 

7= o + 7r= 

' 2^ ^ Nr 2 



(12) 



Let us analyze the simplest case where the masses rui are 
degenerate. We will assume symmetric boundary condi- 
tions for the matter fields. To be more precise, we will 
consider Wi = —1 (and therefore k = Nf) for the path 
drawn by the domain wall. We can then assume that 
all the matter condensates have the same z dependence. 
Then the configuration is described by four real functions 



M^(2) = |M,|p(z)e'"W , 
S{z) = \S,\ Riz)e'l^^'^ . 



(13) 



Notice that we have defined p{z),R{z) in such a way 
that p{±oo) = R{zLoo) = 1. On the other hand, a varies 
from to 2Ti{Nf /Nc - 1) and (3 from to 2TiNf/N^. A 
consistent ansatz under reflection z —> —z is given by: 
p{z) = p{-z), R{z) = R{-z), fi{z) = 2tt^ - (}{-z) and 

a{z) = 2tt{j^ ^ 1) ^ a{—z). Then, we have the following 
boundary conditions at z = 

«(0)=.(f -1), m = . (14) 

Eqs (|l^) imply the following BPS constraint 

Im[e'''WiS,]Vri)] ^ const . (15) 
In particular, at 2 = we have 



Ro [{N, - Nf){lnRo - 1) + Nf Inpo] - A^/Po 

— Nr COS TT 

Nr 



(16) 



where Rq = R{0) and po — p(0). 

The case with Nr = Nc — 1 has already been con- 
sidered for SU{2) [gjS, SU{3) ||l^, and generahzed to 
arbitrary SU{N) in |ll|,|l|l. Since k = Nc — 1, the corre- 
sponding domain walls connect a minimun and its nearest 
neighbour. In these papers it was shown that these do- 
main walls are BPS states only for squark masses lower 
than some critical value, m*, that depends on Nc and 
the Kahler potential. The existence of this bound is re- 
lated to the presence of two different BPS domain wall 
solutions for small enough values of m, which became 
identical at the critical value. 
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FIG. 1. 

p(z) as defined in Eqs (|l^ ) versus z (in units of A~^), for m = 2 
(dotted line), 20 (dash-dotted), 100 (dashed), 200 (solid). 

We have done a similar analysis for other values of Nf , 
using the same Kahler potential , i.e. K. = {SSY^"^ + 
{MMy/^. Here we have worked in detail the case 
Nc = 3, Nf = 1, whereas other cases will be presented 
elsewhere jl^. We have found that the equations can 
be solved for all values of the squark mass, and we have 
checked that the logarithmic branch is never crossed. The 
profiles for p and R are shown in Figs |l|,|| for several val- 
ues of m (given in units of A), focusing on their central 
region. The spatial coordinate z is expressed in units of 
A-i, where A = A(^)^//3^'= is the effective QCD scale 
that arises in the large m limit. 

In our case there is only one BPS solution for every 
value of m. This can be understood analyzing both the 
large and small m limit: 

• When m << A and if 55 is non singular, we can inte- 
grate out the gaugino condensate by imposing dW/dS = 
0. The corresponding Wess-Zumino model describing the 
matter condensate has a BPS state with the required 
boundary conditions for all values of Nf < Nc- In the 
case analyzed by Smilga and Veselov there is yet 

another BPS solution that cannot be described by inte- 
grating S out and that corresponds to S* ~ 0. The exis- 
tence of this domain wall is probably related to the fact 
that the Kahler metric if 55 is singular at this value [ p^ . 
In this limit S ^ 0, po can be derived from Eq. (|l|), and 
it is given by 

Nc f Nf\ 

In our case, since Nf/Nc < 1/2, the resulting value for 
Pq is never positive and therefore we do not find a second 
branch of solutions to the BPS equations. 



FIG. 2. 

R{z) as defined in Eqs (^^ versus z (in units of A~^), for m = 2 
(dotted line), 20 (dash-dotted), 100 (dashed), 200 (solid). The thick 
solid line corresponds to the m ^ oo solution given by Eq. (^o|). 

• Let us turn now to analyze large mass values, m >> 
A. From Fig. ^ we see that there is a well defined gaug- 
ino condensate profile in the m 00 limit. In fact, if we 
assume that this limit exists, the following constraints 
should apply in the asymptotic regions 

p(z)e^"(") = R{z)e"^^'-^ {z « 1/m) , 

(18) 

/9(z)e'"(") = i?(z)e*('^(")-2.) (z » 1/m) . 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 



FIG. 3. 

Plot of the combination ^(/3(z) — 0(2)) versus z (in units of A^-^), 
form = 2 (dotted line), 20 (dash-dotted), 100 (dashed), 200 (solid). 

In Fig. ^ we have drawn the combination ■^{(3{z) — 
a(z)), confirming the previous statement; also a quick 



glance at Figs 1,2 tells us that both p and R follow iden- 
tical paths in the asymptotic regions. 

Therefore, using Eq. in this large m limit we can 
get rid of a and p in the BPS equations. Also the BPS 
constraint involves only the gaugino condensate and can 
be written as 



= const 



(19) 



where I3{z) = (3{z) for z < and /3(z) = P{z) - 2TTNf/Nc 
for z > 0. This constraint allows us to express /? as a 
function of R, and we end up with the following BPS 
equation for R{z) 

d,R{z) = 6iVc(i?(z))'*/3A |cos(7 + (3[R{z)]) lni?(z) 

- sin(7 + /3[i?(z)])/3[i?(z)]} , (20) 

together with the boundary condition at the origin cal- 
culated from Eq. ([l^), i.e., 



i?o(l-lni?o) =cos (^^^^ 



(21) 



As we can see from Fig. 4, this equation has always a so- 
lution with i?o > 1. There is also a solution with i?o < 1 
when 2Nf < Nc- We have seen that only the case i?o < 1 
gives a finite, domain wall like profile. We have also ver- 
ified that the corresponding solution for R{z) is precisely 
the large m limit profile, which is represented by a thick 
solid line in Fig. 2. Therefore the condition ensuring that 
there is just one branch of BPS states at low m values 
also guarantees the existence of the large m limit case. 
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FIG. 4. 

Contour plot of the constraint Eq. (^) in the plane defined by the 
variables _Ro (x-axis) and Nf/Nc (y-axis). 



In summary, it is possible to build BPS domain walls 
in SQCD with 2A''/ < Nc, both in the weak coupling 
(Higgs) regime and in the strong coupling limit, where 
the theory approaches pure supersymmetric gluodynam- 
ics. When the Kahler metric is non singular along the 
different configurations the existence of these solutions 
in the strong coupling regime can be understood by just 
analyzing the superpotential, as in the cases we have just 
shown. 
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